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Candidates for detecting exoplanetary radio emissions 
generated by magnetosphere-ionosphere coupling 
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ABSTRACT 

In this paper we consider the magnetosphere-ionosphere (M-I) coupling at Jupiter- 
like exoplanets with internal plasma sources such as volcanic moons, and we have 
determined the best candidates for detection of these radio emissions by estimating 
the maximum spectral flux density expected from planets orbiting stars within 25 pc 
using data listed in the NASA/IPAC/NExScI Star and Exoplanet Database (NStED). 
In total we identify 91 potential targets, of which 40 already host planets and 51 have 
stellar X-ray luminosity 100 times the solar value. In general, we find that stronger 
planetary field strength, combined with faster rotation rate, higher stellar XUV lumi- 
nosity, and lower stellar wind dynamic pressure results in higher radio power. The top 
two targets for each category are e Eri and HIP 85523, and CPD-28 332 and FF And. 
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1 INTRODUCTION 

The construction of next-generation radio telescopes 
such as LOFAR has in recent years provided impetus 
for considering which might be the best targets for 
detecting exoplanetary auroral radio emissions. This 
attention has traditionally focused on so-called 'hot 
Jupite r'-type exoplanets orbi ti ng cl os e to their loc a l star 



a process known as the Vasy liunas Cyc l e (e.g. Iffilil [l979l; 



ie.E. |rarreii, uescn 6s z^arKai iiyyy 
Zarka et al. 200 ll, 2007; Lazio et all 


; rarrcu et an i/uum; 
20041; Griefimeier et al. 


2004, 2005; Gricfimcicr, Zarka & Sprceuw 2007; Stevens 


20051; iJardine & Cameron! 20081 


iFares et all l20ld: 


Reiners & Christensen 2010l; Vidotto. Jardine & Helling 



'llD . and in such cases the auroral radio emission is 
assumed to be generated by a star-planet interaction, 
mediated either by magnetic reconnection as at the 
Earth or Alfven waves such as at Io. However, iNicholsl 
|2011al ) recently considered the radio emission generated 
by magnetosphere-ionosphere (M-I) coupling at Jupiter-like 
exoplanets with internal plasma sources, and concluded 
that such systems are also able to generate detectable 
emissions. In this process, shown schematically in Figure [TJ 
plasma is generated internally to the magnetosphere from 
sources such as volcanic moons or ionospheric outflow 
and, in a fast-rotating magnetosphere such as Jupiter's, 
becomes centrifugally unstable and diffuses radially away 
from the planet before being lost down the dusk flank 
of the magnetotail via the pinching off of plasmoids, in 
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| Vasvliunadll98 j lHilj|20^ fc Buncel 

boOll; iNichols fc Cowlevl l2003t 12004 . l2005l ; l^cholsl l20nbl) 

Conservation of angular momentum requires that, as the 
plasma diffuses radially outward, its angular velocity drops 
such that a radial gradient of equatorial plasma angular 
velocity is set up, which when mapped along the magnetic 
field to the planet causes a current to flow in the Pedersen 
layer of the ionosphere. This Pedersen current balances 
through the J x B force the drag of the atmospheric 
neutrals on the sub-rotating ionospheric plasma, and this 
torque is transmitted along the field lines to the equatorial 
plane by the sweep-back of the planet's magnetic field out 
of the meridian planes, such that the Pedersen current is 
balanced by an associated radial current in the equatorial 
plane. Current continuity is maintained between these two 
field-perpendicular currents by field-aligned currents, the 
upward component of which, associated with downward- 
precipitating electron s, is on Jupite r associ a ted with the 
main auroral oval dGrodent et all l2003al ; IClarke et all 
2004'; INichols et al. 2009b) and significant components of 
the planet's radio emis sions, i.e. t he b-KOM, HOM and 
non-Io-DAM emissions (|Zarkalll998l ). 

INicholsl | 2011a) showed that the best candidates for 
detection of such internally-generated radio emissions are 
rapidly rotating Jupiter-like exoplanets orbiting stars with 
high X-ray-UV (XUV) luminosity at orbital distances be- 
yond ~1 AU. The XUV luminosity is significant since it 
generates the ionospheric conductivity which allows intense 
M-I coupling currents to flow. It is also worth noting that 
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Figure 1. Sketch of a meridian cross section through a 
Jupiter-like exoplanet's inner and middle magnctosphcrc, show- 
ing the principal physical features involved. The arrowed solid 
lines indicate magnetic field lines, the arrowed dashed lines the 
magnetosphere-ionosphere coupling current system, and the dot- 
ted region the rotating disc of outflowing plasma. From Nichols 
(2011a). 



the variable nature of very active stars tends to reduce the 
frequency of detection of exoplanets around such stars us- 
ing radial velocity and transit methods. An obvious ques- 
tion is then how many internally-generated radio targets 
might one expect to exist, and we thus address this issue in 
this Letter. We determine the best candidates for detection 
of these radio emissions by estimating the maximum spec- 
tral flux density expected from planets orbiting F-M dwarfs 
within 25 pc using data listed in the NASA/IPAC/NExScI 
Star and Exoplanet Database (NStED, now the NASA Exo- 
planet Archive) . We show that a number of systems may be 
detectable if they host massive, fast-rotating planets. The 
results discussed here will be of benefit for the currently- 
underway planning of observations of exoplanetary candi- 
dates using LOFAR, which will have a sensitivity threshold 
of ~1 mjy at 70 MHz and ~0 .14 mjy at 240 MHz for 1 h 
integration (|Farrell et al.|[2004T ). For simplicity, here we use 
1 mjy as the detection threshold, and note that for higher 
frequencies this is thus a conservative value. 



2 ANALYSIS 

We take the X-ray luminosity Lx as a proxy for the 
XUV band as a whole, si nce X-ray and EUV lu minosities 
are broadly correlated l|Hodgkin fc Pvei [l994!), and we 
employ values of Lx determined from data listed in the 
NStED catalogue. Values of Lx are either those as directly 
measured during the ROSAT all-sky survey (RASS), 
including flux detected in both the hard (0.5 -2.0 keV) and 
soft (0.1 - 0.4 keV) passbands of ROSAT l|Hunsch et all 
1999), or calculated from th e chromospheric emissio n 
ratio -Rhk using the relation of ISterzik fc SchmittJ ||l997l) . 
either using the directly-measured value or that computed 
fro m the Moun t Wilso n S-value using the relations given 
by iNoves et al.l (| 19841 ). These relations are given in the 
Supplementary Material (SM). It is worth noting, however, 
that the stars with the highest XUV luminosity generally 
have directly-measured values. The proxies used here are 
typically related to the activity of the star, parameterised 
by the ratio Lx/L^ a \, although here we are simply inter- 
este d in Lx- We do no t e, how ever, that all stars for which 
the ISterzik fc SchmittJ [jl997r i relation is employed have 



Lx/Lboi < —3.7, i.e. within the range of validity of the 
relati on and not within the activity-rotation saturation 
zone (|Giidelll2004l ). 

The value of Lx determined as above for each star 
within 25 pc (within error) listed in the NStED database 
is then used to determine the Pedersen conductance of the 
planet, which comprises a component decreasing with or- 
bital distance as 1/ R or b and a constant component induced 
by auroral precipitation (Nichols 2011a), i.e. 



J xuv 



1/2 



2.6 

Rorb 



+ 1.5 mho 



(1) 



Details of the computation of the M-I coupling cur- 
rents and radio power are given in previous works (e.g. 
iNichols fc Cowlevil2003l ; lNicholsll20lTal ). but briefly, the cur- 
rents arise from an angular velocity gradient in the magneto- 
sphere owing to the centrifugally-driven outflow of plasma, 
described by 



p e d 
2 dp e 



where p e is the distance from magnetic axis, to is the plasma 
angular velocity, £l p is the planetary angular velocity, M is 
the plasma mass outflow rate (taken here to be the canonical 
jovian mass outflow rate of 1000 kg s _1 ), \B ze \ is the magni- 
tude of the north-south magnetic field threading the equato- 
rial plane, and F e is the equatorial value of the poloidal flux 
function related to the magnetic field via B = (l/p)X7F x ip. 
The equatorial field strength \B ze \ is dependent on the plan- 
etary equatorial field strength B eq and the size of the mag- 
netosphere. Here we examin e results of two models for B Rq . 
The first is that employed bv lNicholsl (|2011al ), who took B eq 
to vary with the rotation rate of the planet as tip 1 ' 75 , and 
we consider planets with £l p /£lj — 1 and 3, i.e. represen- 
tative of the angular velocities that might be expected of 
Jupiter-mass planets, thus corresponding to B eq /Bj = 1 
and ~ 2.3, respectively, where Bj is Jupiter's equatorial 
surface field strength. For comparison, we also consider a 
more massive planet with mass 10M j , for w hich we employ 
the model of iReiners fc Christensenl |2010), which is inde- 
pendent of rotation rate for fast rotating bodies, and for 
lOAfj yields B eq /Bj — 17. These magnetic field strengths 
correspond to polar electron cyclotron frequencies, and thus 
radio emission bandwidth, of ~24, ^54, and ~406 MHz, 
respectively. The size of the magnetosphere is governed by 
pressure balance between magnetospheric magnetic field and 
plasma pressures on the one hand and stellar wind dynamic 
pressure on the o ther. Along with iNichoTsl (|2011al ) we em- 
ploy the empirical iHuddleston etaL ( 19981 ) relation for sub- 
solar magnetopause distance versus dynamic p ressure p S w, 
a qua ntity related to the stellar mass loss rate. IWood et alj 
(2005) have discussed a relation between stellar activity and 
mass loss rate per unit surface area. However, whether this 
relation holds for very active stars remains uncertain, and 
they showed that stars more active than the Sun exhi bit val- 
ues ra nging between ~0. 01-100 times the solar value. ICohenl 
l|201ll ) showed that for the Sun no relationship between mass 
loss rate and X-ray luminosity is apparent, and attributed 
this to the fact that the solar mass loss is dominated by the 
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fast solar wind originating from open solar magnetic flux, 
while solar activity is associated with closed flux. In the light 
of this uncertainty, and in the lack of in situ measurements of 
stellar wind densities and velocities we simply consider two 
stellar wind dynamic pressures, p sw q and lOOp^Q. Once 
the plasma angular velocity profile has been obtained from 
equation [2] the density jiu of the upward-directed auroral 
field-aligned current responsible for the radio emission is 
given by 



d 



p e B ze dp e 



(3) 



This field-aligned current must in general be carried by a 
field-aligned voltage, which yields a precipitating electron 
energy flux given by 



2 



J\\i 

j\\io 



(4) 



where Ef a and are the energy flux and field-aligned 
current that can be carried by unaccelerated precipitating 
electrons. The radio power P r is then found by integrating 
this precipitating energy flux over the hemisphere and 
converting to radio power assuming a 1 per cent generation 
efficiency for the electron cyclotron maser instability. 
The maximum power was determined over the range 
1 < (Rorb/AU) < 500. The outer limit is somewhat arbi- 
trarily large, but we note that exoplanets with semi-major 
axes of up to several thousand AU have b een detected 
by direct imaging (e.g. iBurgasser et al.l 120101 ). Finally, the 
spectral flux density T is determined assuming the em ission 
is beamed into 1.6 sr ijZarka. Cecconi fe Kurthll2004h with 
bandwidth equal to the polar electron cyclotron frequency. 

Computed powers versus orbital distance for the 
different cases discussed above are shown in the SM, but 
in Figure [2] we show the maximum radio powers computed 
by the model versus Lx for the different cases. For each 
case the results employing both terms in equation Q] are 
shown, along with those for only the stellar XUV term for 
comparison. First, comparing panels it is apparent that 
higher magnetic field strengths and faster rotation lead 
to hig her radio powers. It is also clear t hat, as discussed 
e.g. bv lCowlev. Nichols fc Andrews! |2007l ). higher dynamic 
pressure results in lower auroral current intensities, and 
thus lower radio powers, owing to the lower magnetosphere 
size. This effect also results in increased power with orbital 
distance, an effect which is, however, offset by the simul- 
taneous decreasing of the conductance induced by stellar 
XUV flux. Considering then solely the effect due to XUV 
flux, the power reaches a maximum at an orbital distance 
that is larger for increased Lx, resulting in increased 
maximum radio power as shown by the dashed and dotted 
lines in Figure [2] Including the conductance contribution 
from auroral precipitation results in the power asymptoting 
to a finite value at R or b = oo instead of decreasing to zero, 
resulting in the flattening of the solid and dot-dashed lines 
in Figure [2] for low XUV luminosity. For planets with larger 
magnetic fields, the switch between regimes dominated by 
XUV- and precipitation-induced conductance occurs at 
increased Lx, such that the power from planets with large 
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Figure 2. Figure showing computed maximum radio powers P r 
versus Lx ■ The different cases discussed in the text are shown by 
the labels above each panel and in panel (a). In each panel values 
of Lx for stars with planets are shown by the pluses and values 
for stars with Lx > 100 are shown by the crosses. 



field strengths are essentially independent of Lx- In these 
cases the spectral flux density is instead strongly dependent 
on distance from the Earth to the star. 

Two categories of target were considered in the anal- 
ysis. The first is all F-M dwarfs within 25 pc (within er- 
ror) which are known to host planets and have directly- or 
indirectly- measured values of Lx, and the second is those 
with Lx greater than 100 times the solar X-ray luminosity 



irrespective of whether the star is known to host plan- 



ets or not. We take L x q to be th e mean of the 0.1-2.4 keV 
full so lar cycle range measured by I Judge. Solomon fc Avresl 
i|2003h . i.e. 10 20 ' 35 W The threshold of 100L xQ was used 
since in the model of iNicholsl i|2011al ) stars with such high 
X-ray luminosities are required to host planets detectable 
from beyond 1 pc with the jovian rotation rate. 



3 RESULTS 

A total of 40 F-M dwarfs known to host planets have 
directly- or indirectly measured values of Lx, while 51 have 
Lx > 1OOL X 0. A complete listing of the results of this 
study is available in the SM, but in Tables [1] and [2] we show 
the top 10 results for stars known to host planets, and stars 
with Lx > 100L x q, respectively. It is first apparent from 
columns 6 and 8 of Table[T]that no planets of Jupiter's mass 
which rotate at the jovian angular velocity and orbit stars al- 
ready known to host planets would be expected to be viable 
targets. The situation is improved somewhat for the faster 
rotating planets considered in columns 7 and 9, with 15 stars 
potential targets for detection of a planet with (Q p /Qj) = 3, 
reducing to 10 for 100p s „Q, although we note that most of 
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Table 1. Table showing the top ten potential targets for observing internally-generated exoplanetary radio emissions amongst stars known 
already to host planets. The major columns are the SIMBAD default name, the spectral type, the distance s in parsecs, declination in 
degrees, the star X-ray luminosity Lx , and estimates of the maximum radio spectral flux density T . The Lx column indicates whether the 
value is determined from direct ROSAT measurements (R) or indirectly from the Mount Wilson S-value (S) or chromospheric emission 
ratio (H). For each target, eight estimates of the maximum radio spectral flux density T are shown. Four of these employ the 

planetary magnetic field strengths used by Nichols (2011), and four employ the Reiners et al. (2010) algorithm for a planet with mass 
M p = lOMj. Each of these groups are divided into results employing solar and 100 X solar dynamic pressure values, and these are 
further divided into results employing two values of the planetary angular velocity normalised to Jupiter's rotation rate f2 p /f2j. 



T I mjy 

B eq Nichols (2011) B eq Reiners et al. (2010), 10Mj 
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V* eps Eri 


K2V 


3.2 


-9.5 


21.32 (R) 


0.61 


18.21 


0.49 


11.63 


31.43 


282.87 


5.10 


45.93 


HIP 85523 


M2+V 


4.5 


-46.9 


20.53 (R) 


0.30 


9.04 


0.24 


5.78 


15.62 


140.56 


2.55 


22.95 


V* IL Aqr 


M3.5V 


4.7 


-14.3 


19.49 (R) 


0.28 


8.38 


0.23 


5.37 


14.49 


130.42 


2.37 


21.35 


LHS 3685 


M2/3V 


4.9 


-49.0 


19.77 (R) 


0.25 


7.59 


0.21 


4.86 


13.13 


118.19 


2.15 


19.33 


LHS 349 


G5V 


8.5 


-18.3 


19.87 (R) 


0.08 


2.55 


0.07 


1.64 


4.42 


39.75 


0.72 


6.50 


HR 7722 


K3V 


8.8 


-27.0 


20.23 (R) 


0.08 


2.39 


0.06 


1.53 


4.13 


37.15 


0.67 


6.07 


LHS 311 


G3/5V 


9.2 


-40.5 


19.78 (R) 


0.07 


2.17 


0.06 


1.39 


3.75 


33.78 


0.61 


5.53 


LHS 310 


M3V 


10.2 


26.7 


19.84 (R) 


0.06 


1.77 


0.05 


1.13 


3.06 


27.56 


0.50 


4.51 


LHS 3257 


M1V 


10.3 


25.7 


19.25 (S) 


0.06 


1.74 


0.05 


1.11 


3.00 


27.03 


0.49 


4.43 


HD 13445 


K1V 


10.9 


-50.8 


20.61 (R) 


0.05 


1.57 


0.04 


1.00 


2.71 


24.36 


0.44 


3.98 



these stars are in the southern hemisphere. For the WMj 
planets in columns 10-13, 21 would be detectable with a so- 
lar wind dynamic pressure and jovian rotation rate, reduc- 
ing to 4 for 100p su ,Q, while for (Q p /Qj) = 3 all 40 would 
be detectable with solar wind dynamic pressure, reducing to 
34 for 100p sm Q . The brightest target in the northern hemi- 
sphere is LHS 310. Considering now the targets listed in 
Table [2] we first note that the nearest star is at ~8.7 pc, 
i.e. further than the top 5 targets in Table [T] However, for 
the regimes in which the XUV-induced conductance is not 
negligible (i.e. for columns 6-9, corresponding to panels (a) 
and (b) in Figure [2} the higher XUV luminosity results in 
higher radio powers, and thus higher spectral flux density for 
a given distance from Earth than in Table [1] but otherwise, 
the pattern is similar. Thus, again, from columns 6 and 8, 
it is apparent that no planets with jovian mass and rotation 
rate would be detectable, but columns 7 and 9 indicate that 
12 fast-rotating planets would be detectable with solar wind 
dynamic pressure, reducing to 4 for 100p su) Q. For the lOMj 
planets in columns 10-13, 20 would be detectable with a so- 
lar wind dynamic pressure and jovian rotation rate, reducing 
to none for W0p sm Q, while for (fi p /S7j) = 3 all 51 would be 
detectable with solar wind dynamic pressure, reducing to 39 
for lOOp s „ . 



4 SUMMARY AND DISCUSSION 

In this paper we have considered which might be the 
best targets for the discovery of internally-generated ex- 
oplanetary radio emissions. We have employed the stellar 
X-ray luminosity data available in the NStED database 
(now the NASA Exoplanetary Archive) to determine the 
maximum radio power available from M-I coupling due 
to outflowing internally-generated plasma for 2 values of 
each of the planetary rotation rate, mass, and stellar wind 
dynamic pressure. We have considered two categories of 



potential targets: all F-M dwarfs within 25 pc (within 
error) which are known to host planets and have directly- 
or indirectly-measured values of Lx, and all those with 
Lx > 100L x q, irrespective of whether they are known 
to host planets or not. We have shown that up to 40 and 
51 potential targets exist, respectively, for each of these 
categories with the actual number depending on the system 
parameters. In general, stronger planetary field strength, 
combined with faster rotation rate, higher stellar XUV 
luminosity, and lower stellar wind dynamic pressure results 
in higher radio power. The top two targets for each category 
are e Eri and HIP 85523, and CPD-28 332 and FF And. All 
these are in the southern hemisphere. The top two northern 
hemisphere targets in each category are LHS 310 and LHS 
3257, and V834 Tau and OU Gem. 

It is worth mentioning that this model requires a num- 
ber of elements in place to produce a detectable system, such 
that, while we have highlighted eight targets above, a wider 
survey of reasonable targets would increase the chance of 
discovering a system with all element s in place. It is also 
worth re-iterating that iNicholsl l|2Qllal ) highlighted a num- 
ber of areas in which the model could be developed. For 
example, the model does not take into account the stretch- 
ing of the magnetic field due to the ce ntrifuga l force and 
hot plasma pressure, shown recently by INicholsl (|2011bl ) to 
have a significant effect on the magnitude of the currents, 
and we have not considered at all the stellar wind interac- 
tion, which could be asso ciated with signif i cant emissions 
in Jupiter's polar regions llWaite et al. I l200ll ; I Grodent et al.l 
l2003bl ; lNichols et al.ll2009al lbl) . Therefore, it should be noted 
that the list of the best targets for observation may be rea- 
sonably expected to evolve with future developments of the 
model. Finally, while we have not discriminated by multi- 
plicity in th i s stud y, we note that recent observations by 
iDovle et al.l (|201ll ) have shown that exoplanets can exist 
around binary star systems. 
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Table 2. As for Table [T] but for any star with L x > 100L X& . 
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CPD-28 332 


F9-V 


8.7 


-28.2 


22.66 (R) 


0.19 


3.31 


0.07 


1.71 


4.64 


41.77 


0.73 


6.54 


V* FF And 


MOVEP 


8.6 


-20.6 


22.38 (R) 


0.15 


2.84 


0.07 


1.69 


4.59 


41.34 


0.73 


6.56 


V* AT Mic 


M1VE 


9.9 


-31.3 


22.74 (R) 


0.15 


2.55 


0.06 


1.32 


3.57 


32.14 


0.56 


5.03 


V* AK Pic 


M4.0V 


10.2 


-32.4 


22.55 (R) 


0.13 


2.25 


0.05 


1.23 


3.34 


30.03 


0.52 


4.72 


V* BY Dra 


MOVP 


11.5 


-43.8 


22.60 (R) 


0.10 


1.78 


0.04 


0.97 


2.64 


23.72 


0.41 


3.73 


V* V834 Tau 


GOV 


12.8 


47.7 


22.70 (R) 


0.09 


1.55 


0.03 


0.80 


2.17 


19.51 


0.34 


3.05 


HIP 61941B 


F1V+F0 


11.8 


-1.4 


22.42 (R) 


0.08 


1.51 


0.04 


0.90 


2.45 


22.05 


0.39 


3.50 


V* OU Gem 


GOV 


21.7 


33.9 


23.66 (R) 


0.10 


1.36 


0.01 


0.33 


0.92 


8.24 


0.13 


1.17 


V* V775 Her 


K4Vke 


13.5 


20.9 


22.57 (R) 


0.07 


1.29 


0.03 


0.71 


1.92 


17.27 


0.30 


2.72 


* tet Boo 


F8V 


14.6 


51.9 


22.63 (R) 


0.06 


1.11 


0.03 


0.61 


1.64 


14.80 


0.26 


2.33 
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